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ABSTRACT 


The  radiation  impedance  of  a  finite  array  of  finite  cylinders  is 
determined  by  a  boundary  value  problem  technique.  The  final  equations 
for  the  radiation  impedance  are  stated  without  restriction  on  oylinder 
length  or  radius,  oylindor  spacing,  or  on  the  amplitude  and  phase  of 
the  radial  velooity  of  the  cylinders.  The  derivation  is  made  for  the 
case  of  each  cylinder  in  tho  array  having  a  uniform  radial  velocity  over 
a  soctor  of  the  oylinder,  but  tho  method  of  solution  is  not  limited  to 
such  a  velocity  distribution.  The  results  of  a  numerical  calculation  of 
tho  effect  of  electrical  beam  steering  on  the  radiation  impedance  of  a 
specific  array  are  shown  to  illustrate  such  effect  and  the  uso  of  the 
derived  equation  for  radiation  impedance. 

The  writer  wishes  to  express  his  appreciation  for  the  assistance  and 
onoouragoment  given  him  by  I.Ir.  John  l'liokman  of  the  U.  S.  Ilavy  Electronics 
Laboratory  and  Professor  II.  liedwin  of  the  U,  S.  Naval  Postgraduate  School 
in  this  investigation. 
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1#  Introduction 


The  radiation  impedance  of  a  finite  array  of  finite  cylinders  in  an 
infinitely  long,  rigid,  cylindrical  baffle  has  been  investigated  by  D.  H# 
Robey# ^  Robey  considered  the  problem  of  uniformly  pulsating  cylinders 
and  his  final  conclusions  are  based  on  a  cylindrical  array  of  radius 
negligible  in  comparison  to  the  wavelength  of  the  radiated  sound#  It  is 
the  intent  of  this  paper  to  investigate  the  radiation  impedance  of  such  a 
finite  array  of  finite  cylinders  with  fewer  restrictions  imposed  on  the 
radial  velocity  distribution  and  size  of  the  array#  It  is  also  the  intent 
of  this  paper  to  apply  the  results  of  such  an  investigation  to  the  study 
of  the  effects  of  electrical  beam  steering  on  the  radiation  impedance  of 
a  finite  array  of  finite  cylinders# 

The  radial  velocity  distribution  on  each  pulsating  cylinder  will  be 
specified  as  a  boundary  condition  to  the  wave  equation  which  is  to  be 
solved#  However  this  is  not  enough  to  completely  specify  a  unique  solu¬ 
tion  to  the  wave  equation#  Therefore  it  is  proposed  that  as  the  radial 
distance  from  the  cylindrical  radiator  increases,  vanishing  boundary 
values  for  pressure  and  velocity  aro  required#  This  condition  again 
proves  too  weak  to  determine  a  unique  solution#  Therefore  a  more  stringent 
boundary  condition  in  addition  to  the  condition  of  vanishing  boundary 
values  is  required,  ie.  nThe  Sommerfeld  Condition  of  Radiation#*1  That  is, 
"The  energy  which  is  radiated  from  the  souroes  must  scatter  to  infinity; 

^D#  H.  Robey,  On  the  Radiation  Impedanoe  of  an  Array  of  Finite  Cylinders, 
J#  Acoust#  Soo#  Am#, 27,  pp706-711,  1955 

^A#  Sommerfeld,  Partial  Differential  Equations  in  Physics,  Academic 
Press  Inc.,  1949,  p  189 


1 


no  energy  may  be  radiated  from  infinity  into  tho  prescribed  singularities 
of  the  field  (piano  waves  are  excluded  since  for  them  even  the  condition 
u  -  0  fails  to  hold  at  infinity).11 

The  wave  equation  which  is  to  be  satisfied  is  stated  in  cylindrical 
coordinates : 

*  2  r  ^  T  2  2*2  2  ,  .  2 
Or  r  or  roe  0  z  cat 

The  scalar  velocity  potential  function  is  understood  to  be  ^(r,e,z,t). 

The  z  axis  of  the  cylindrical  system  corresponds  to  tho  longitudinal  axis 

of  the  cylindrical  array  (see  figure  l). 

The  driving  forces  which  are  to  be  imposed  on  the  elements  of  the 

cylindrical  array  will  vary  harmonically  with  time.  Therefore  for  the 

results  of  this  paper  to  conform  to  standard  electrical  engineering 

practice  of  assigning  positive  reactance  to  inductance,  a  time  dependence 

of  the  form  will  be  assumed.  That  is  j^(r,  e,z,t)  -  0(r,e,z) 

After  substitution  for  ^(r,ez,t),  tho  wave  equation  is  reduced  to  the 

following: 

4-1 .  4  £&  4  —1 _ 4  £$  z  44L.  <a 

>  2  r  or  2  -.2  *2  2 

dr  r  de  dz  c 

A  formal  solution  to  the  above  equation  is  discussed  in  Appendix  1. 

The  radiation  impedance  of  an  array  of  finite  cylinders  will  now  be 

determined. 
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20  Derivation 


The  first  array  to  be  investigated  is  the  least  restricted  in 
boundary  conditions  of  the  arrays  to  be  consideredo  The  array  is  to  be 
constructed  as  follows  (soe  figure  l): 

1.  There  are  m  active  cylinders  of  length  =  2^*,  radius  =  a, 
and  center  to  center  spacing  =  2b* 

2 •  The  cylinders  are  mounted  coaxially  in  an  infinitely  long, 
rigid,  cylindrical  baffle  of  radius  =  a# 

The  pulsation  of  any  active  cylinder  is  to  be  described  as  follows: 

1*  On  the  given  cylinder,  the  radial  velocity  is  constant 
over  the  length  of  the  cylinder  for  a  given  aximuth  angle* 

2.  On  all  cylinders,  the  radial  velocity  is  uniform  between  azimuth 
coordinates  -  <X  to  */  oC  and  zero  elsewhere* 

Therefore  the  radial  velocity  of  the  i^*1  cylinder  may  be  described  by  the 

following  Fourier  Series: 

oo 


h  - 


(2o0l) 


Where  A  s  sn  sin  nc< 
n 

>  en 

=  1» 

o 

D 

£ 

ncc 

6n 

=  2, 

n  “  1,2,3, •••• 

and  -  the  phase  of  the  pulsations. 

The  scalar  velocity  potential  0(r, e, z,t )  must  satisfy  the  following 
wave  equation: 

2 

V2  0  =  h  ^  f  (2.02) 

0  TT* 
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The'scalar  velocity  potential  is  assumed  to  vanish  as  r-^oo  and  the 
Sommerfeld  Radiation  Condition  is  assumed  to  bo  valid  vrLth  respect  to  the 
radial  coordinate* 

The  time  dependence  of  the  scalar  velocity  potential  is  taken  to  be 
=  tf{r,  e,z)eJu,t  as  discussed  in  Section  I*  The  scalar  velocity 
potential  $(r9e,z)  must  then  satisfy  the  Helmholtz  V/avo  Equation 
(V^/k^)  fS  -  0*  (2*02)  An  eigenfunction  which  is  a  solution  to  the 

Helmholtz  Wave  Equation,  satisfies  the  vanishing  boundary  conditions  and 
the  Sommerfeld  Radiation  Condition,  is  developed  in  Appendix  1* 


This  eigenfunction  is  stated: 

d  (  \  -  a  sin(ne)  =j2TT  J  j  (2) 

K?  ~  V?  cos(ne) 


v  Vd 


Where 


y«2  a  k2  -  (2 n  f  )z 


(2.03) 

(2.04) 


/  =  l|k2  -  (2  Trf  )2’ 


(2  7rf)z<  k2 
(2  ff  f)z  >  k2 


^  =  -j  filr  r  )2  -  k2" 

It  is  now  proposed  to  determine  a  $  ( r,e,z ),  constructed  from  a  sum 
of  the  eigenfunctions  (r,e,z),  which  will  satisfy  the  boundary  condi- 

tions  at  the  surface  of  the  cylindrical  array*  The  technique  to  be 
employed  in  determining  the  complete  scalar  velocity  function  is  a 
logical  development  from  the  technique  of  L*  V*  King  used  in  the  solution 
of  the  radiation  impedance  of  a  piston  in  an  infinite  rigid  baffle*  The 
constant  A^^  transforms  into  *n(  f  )  as  the  eigenvalue  f  becomes  a 
continuous  variable* 


The  sum  of  eigenfunctions  is  taken  to  be: 
^o© 
oo 

0(r,e,z)  a 


n°0 


^  cos(ne)  Hn(2)(/r)  ejZ7rT  2  df 


05) 


-oo 


°L.  V.  Kihg,  On  the  Acoustic  Radiation  Field  of  the  Pie j:6-Electric 
Oscillator  and  the  Effect  of  Viscosity  on  Transmission,  Can.Journ,  Res.,  11, 
pp  135-155,  July-Dee .1934 
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The  velocity  in  tho  radial  direction  iss 
vr(r,ez)  =  -  -fr  0  (r,e,z) 


oo 


n=0 


(2)  (2) 

fn(J>)  Wl  (A 7  ^rVr>  cos(ne)  .W* 

(2.06) 


The  radial  velocity  distribution  of  the  finite  array  will  now  be 
investigated#  From  the  boundary  conditions  as  defined,  it  is  possible  to 
write  an  expression  for  the  radial  velocity  amplitude  at  the  surface  of  the 
array: 


vr(a,e,z) 


oo 


n=0 


g(z)  An  cos (n® ) 


(2.07) 


The  funotion  g(z)  is  equal  to  zero  on  the  rigid  baffle  and  on  the  i^*1 

cylinder  is  equal  to  e^i.  If  the  two  expressions  for  the  radial 

velocity,  (2.06)  and  (2.07),  are  compared;  the  equality  of  the  two 

quantities  is  possible  if  they  are  equal  term  by  term,  that  is: 

(2)  (2) 

eUjAn.  (Ar)  ^ 

(2.00) 


-OO 


Define  a  new  quantity. 


(2)  (2) 

^■/l(|Aa)  ~  ^n-1  (,Aa ) 

2 


A 


(2.09) 
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Then 


-\00 


g(z) 


g(  Y  )  o  -/jZfry  z  d  y 


(2.10) 


-CD 


It  can  new  be  seen  that  G( f  )  and  g(z)  are  a  Fourier  Transform  Pair. 
Therefore : 


G(  f )  «  \  g(z)  e"J27r^  z  dz 


(2.11) 


-oo 


'b-yS 

(0)  Q-iZir$  z  dz  -f  \  V1  eJ^l  ^  z  dz 

-oo  ^  "b-fZ 

(2i-l)b-/,S 

{  Vi  e^i  e-J'27rlf  z  dz  / 

(2i-l)b-/Q^ 

p  (2m-l  )b-//3 


^ . / 


N  00 


Vmej^m  e  -J2rry  2  d2  ^ 


(0)  e"j25r^z 


(2m-l  )b-yS 


(2m-l)b-/^ 


G(J)  ;  vx  1.-32"- ]T  1>  -f 


V 

V  Vm  e-J2Tr^  (2m-l)b 

vT 


m 


G(J)  = 


i=l 


V.  e“027r^  (2i-l)b 


sin  2JT%£> 

^jr 


dz 


6 


(2.12) 


It  is  now  established,  after  substitution: 

m 


?nW  = 


2An  sin(Ziffp)  V±  e^1  e"j2  ^ ?  (2i“l)b 

i=l 


- ru) — — no — ~ 

LVl^'a)  "  Hn-l^a)J 


(2.13) 


and 


/( r,e,z)s 


oo 


n=0 


oo 


X 


V* 

L-oo 


(2) 

2An  cos(ne)  siu(2^  )Hn  i^x) 

rt2i  i2)  v 

/“MT  |Vi^a>  -  Vi/a>] 


m 

Vle^i  e  "J2 (2i-l)b  ej2/Tjf»  z 

its  1 

(2.14) 

To  obtain  the  average  pressure,  one  uses  the  force  equation 

vp  /far  :  0 

Sinoe  V  ■  -  XJ  ft  and  for  the  condition  of  time  dependence  e^1^,  the 
force  equation  takes  the  form: 

p(r,e,z)  «  jw^(r,e,z) 

The  average  pressure  amplitude  on  the  active  area  of  the  hbb  cylinder 
defined  as  p^. 

v/°c  (2h-l)b+y3 

0(a,e,z)  de  dz 


Ph=  J 


•  Co 


4°tJ3 


-cx.  (2h-l  )b-yS 


(2.15) 


is 
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The  averaging  process  is  carried  out  and  the  expression  for  the 


average  pressure  is,  upon  substitution  for  An: 

oo 


oo 

Ph=jw^°  ^ 


n»0 


^2 


sm 


(n°0 


11  (n4 


2  H 

-n . 7  -  ■ 

^2  _  (2) 


(2  7 rtP  )  /[Vl>a)'Hn.l^5 


-OO 


oo 


V.  (i“h)b 


x  s  e^'  •■•e 


i-1 


dl 


(2.16) 

The  average  radiation  impedance  per  unit  area  of  the  active  area 


of  tho  hb^  cylinder  will 

be  defined  as 

<z>  h* 

<Z>h  = 

ik 

V  OYh 

h  e 

Further  define  Q~  a 

IT/4 

<r  z 

^i  -  (*y  )2’ 

(Vf  )2  <  1 

<r  z 

j-  )2  ■ 

•  1  (Vf  )2  >  1 

(2.17) 


(2.18) 

\ 

Then  the  radiation  impedance  per  unit  area  of  the  active  area  of  the 
th 


h  cylinder  may  be  written: 


oo 


<2>h=i4fi»0^ 


n=0 


sin  (n<x) 
'n  3 

(n  «)2 


P  ( 2 ) 

sin  (2  (fl~  ka) 


-OO 


8 


X 


m 


V, 


e-34>f  (i-h)b  djo 


(2.19) 


The  radiation  impedance  of  an  array  of  finite  cylinders  which  pulsate 
uniformly  with  arbitrary  amplitude  and  phase  over  a  sector  of  each  cylin¬ 
der  has  been  determined*  Evaluation  of  equation  (2*19)  involves  numerical 
integration  of  an  infinite  integral*  The  integrand  of  (2*19)  must  be 

carefully  examined  in  the  region  of  -i-*  A  further  handicap  to 

A* 

evaluation  of  (2*19)  is  the  fact  that  second  order  Hankel  functions  of 
complex  arguments  must  be  evaluated*  These  qualities  of  (2*19)  are 
especially  undesirable  if  tho  radiation  impedance  is  to  be  evaluated  by 
digital  computer  techniques*  Equation  (2*19)  will  now  be  reduced  to  a 
form  suitable  for  numerical  integration* 

The  following  relation  is  stated  and  used  in  deriving  a  new  form 
of  equation  (2.19)*4 

Kv( jx)  a  -  i-ire"ji7r  (v_1)  hv(2)(x) 

1 


oo 


n«0 


sm 


H- 


(2) 


■n 


( CT^ka  ) 


!  ,  r  (2)  (2) 


m 


i=l 


vi  „ j (V'i-lh) 


cos  £  4  It  ^  b( i-h )J 


McLachlan,  Bessel  Functions  for  Engineers,  Oxford  at  the  Clarendon 
Press,  1955,  p.  204. 
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oo 


oo 


■ti*  ¥■  2  Sn 

n=0 


i 

1_° 

-K 


sin^(2  Tk j* /3  )(-j)  Kn(.i<Tka)  x 

(2^(f/3)2  [x^  (j<Tka)-/Kn.1(jrka)] 


m 


X  /  -h-  eJ  )  oos  |-47r^  b(i-h)]  d  f 


i=l 


(2,20) 


The  integrals  will  now  be  considered  separately.  Let 


1 


II  = 


sin2  (2^^)  42^(fka) 


0 


o  r(2^   <2) 

(2 irjp)2  <r  [Hn^i^ka)  -  H^^ka)] 


oos  [4^^  (i-h)bj  d 


(2.21) 


Let  s  sin  £ 

then  when  J:  0  ;  j  =  0 

X=f;  f-iw 


«"  =  -  (  *■?  )2'  =  oos  f 

dj  a  COS  ^ 


A. 


! 


Therefore ; 


sin2(k^3  sin  |  ) 

(k^sinj  )2  Hn^(ka  COS  /)  -  Il<5}(ka  COS  /) 


^n^2^(ka  cos J  ) 


0 


cos  2kb ( i-h ) s in jj d  ^ 

~ 


(2.22) 
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Consider : 


oo 


.  2 


I2  = 


1 

x 


Let 


in  (Z'y.f/S)  (-.1)  Kn(,i<rka) _ 

(Z7rf?  )V  [i^  (jr  ka)  T^n-lQ'*'  ka)j 


cos  [4irJ(i-h)bJ 


X'f  B  O  SC  ^ 


Then  -when  -  1_  ;  ^  a  -g-H' 

X 

J  =  oo  ;  /  5  0 


<r=  -j  Jujf  )2  -  T1  =  -j  cot/ 


dj  =  -esc/  cot  f 


.<  <LL 

A 


sin^(k g  esc  /  ) 


Kn  (ka  c  ot  /  ) 


(kyO  )2  cso/  (ka  cot  f  )  Kn-1(ka  cot/  ) 


(2.23) 


cos  ^  2kb  (i-h)c  sc/J  d^ 


(2.24) 

The  expressions  (2.22)  and  (2.24)  are  substituted  in  (2.20)  and  the 
resulting  equation  for  radiation  impedance  per  unit  area  z  ^ takes  the 
more  tractable  form: 


-!r 


oo 


<z>  h  3  J®(°oG  ■p 


n=0 


sin^(noc  )  \  sin^(k/sin/)  n 


H  ^^(ka  cos  £  ) 
n  0 


(n«,)‘ 


(k/,SinS>)^-oSp-.l<!/ka  oosj) 
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m 


cos  [2kb(i-h)sinj]dj 


i=l 


oo 


/  jQ(c 


n=0 


h 

V 

2 


c K,i  sin2  (not.) 


sin^ (k/3  esc  ^  )  Kn(ka  cot  ^  ) _ 

(k^6  )  esc  f  K^-^ka  oot£)  /Kn-^(ka  cotj) 


m 


i=l 


(2*25) 


The  evaluation  of  <z>  h  ^as  b0en  reducod  to  that  of  evaluating  finite 
integrals  and  functions  of  positive  real  variables  for  which  mathematical 
tables  are  readily  available#  A  further  manipulation  of  (2*25)  will  be 
made  in  order  to  separate  the  real  and  imaginary  components  of  < z>  ^# 

Let 


Bn(ka  cos  ^  )  e“j^n(ka  cos  £  ) 


IIn(2)(ka  cos  /  ) 


II 


(2)  ,  (2) 

y^(ka  cos  ?  )  -  Hn_^(ka  cos^  ) 

(2.26) 


Then  the  expression  for  < z>  ^  transforms  into  the  following  expression: 
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<2V  ^8(6°  Jx 


oo 

X 


n=0 


oo 


o  sin2(n  ) 
'n  ,  *2 

(not) 


■  *in  (fr/2i.?*n  O  BQ(ka  cos  |  )^>  -i-  sinP^ka  cos£)  -  (^i“f^)  cos^2kb(i-h)sin^  j 

J  <k/,«lnP  ^vi>  L° 

0 


-8^°  W 


n=0 


oo 


i-Tr 


„  sin2(nc<)  \  sin2(kp  cso  ^  ) 

Cn  2  \  ^2  ^ 

(nc<) 


KR(ka  cot  £ ) 


^  (k/s)  esc  ^  Kn7Zi  (ka  cot^  )/Kn>i(ka  cot!^ 

I 


m 

-i-  sin(^-^)cos  j^2kb(i-h)csc£  j  df 

>4V» 


8<^  ff  2 

n=0 


oo 

/j  8fo°  ^ 

n=0 


hrr 


.  sin2(noc  )  \  sin^' 

n  .  ,2  y 

(n<*  )  \  (k/5  sin  j  ) 


I  IK 


«in2(no<)  \  sin2 

€  - “ 


m 


^  ^.AiELj ,  1  B  (ka  cos  <:  )  ^>  -i-  cos  (ka  cos  ^  )  -  (^-^)cos  T2kb(  i-h)  sib  ^  "j  d£ 
(k/  «inj)2  n  "  ^  Vh  n  - 


(not)2 


\  — ~ - iJJ -  oos(^.-^)cos  |2kb(i-h)csc/  ]  d£ 

\  CSCS  ^-/l  (ka  cot  j  )-/Kn_i(ka  cot^)^  h 

oJ 


13 


The  radiation  impedance  per  unit  area  of  an  array  of  m  finite 
cylinders  which  pulsate  uniformly  over  a  given  sector  of  each  cylinder 
with  phase  and  amplitude  arbitrary  but  constant  for  a  given  cylinder  is 
expressed  by  equation  (2.27).  There  are  no  restrictions  stated  on  the 
wavelength  of  the  radiated  sound,  on  the  cylinder  length  or  radius,  or 
on  the  cylinder  spacing#  Although  equation  (2.27)  is  for  a  system  with 
uniform  velocity  and  phase  across  each  radiating  element,  the  method 
of  solution  is  not  limited  to  such  a  velocity  distribution.  Any 
velocity  distribution  in  azimuth  which  can  be  expressed  in  a  Fourier 
Series  is  possible  of  solution  using  such  technique.  Further  any 
velocity  distribution  in  the  z  coordinate  for  which  a  Fourier  Transform 
exists  is  also  capable  of  yielding  a  solution.  The  distribution  in  the 
z  coordinate  need  not  be  tho  same  on  each  cylinder  for  a  solution  could 
be  constructed  piece  by  piece  for  such  a  system. 
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3.  Spooial  Cases 


The  validit/  of  equation  (2, 27)  will  be  considered  by  investigating 
two  special  cases  of  the  finite  array  for  which  solutions  already  exist. 
The  first  array  to  be  considered  is  the  single  uniformly  pulsating 
cylinder  in  an  infinitely  long,  rigid  baffle*  Let  the  array  be  con¬ 
structed  as  follows: 

1*  There  is  one  active  cylinder  of  length  -  2  f>  and  radius  =  a. 

2.  The  cylinder  is  mounted  coaxially  in  an  infinitoly  long,  rigid, 
cylindrical  baffle  of  radius  =  a. 

3.  The  cylinder  pulsates  uniformly  with  radial  velocity  = 


V0 


The  Fourier  Series  equivalent  of  the  third  condition  as  stated  above 
is  expressable  in  one  term.  That  is,  on  the  active  cylinder  tho  radial 
velocity  may  be  expressed  as  one  term  of  the  Fourier  Series  (2.0l). 


With  reference  to  equation  (2.27),  tho  above  statement  or  condition 
reduces  the  summation  with  respect  to  n  to  one  term,  namely  n  =  0.  The 
condition  of  one  cylinder  reduces  the  summation  with  respect  to  i  to  one 
term,  namely  i  =  0©  Let  the  radiation  impedance  per  unit  area  of  the 
uniformly  pulsating  cylinder  be  <  z>  •  Then  from  equation  (2.27), 


0 
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JT 


\ 


•2 


■/ ^ 


0 


(3.01) 


Equation  (3.01)  is  written  with  substitutions  according  to  (2.26): 


T 


0 


(3.02) 


0 

Equation  (3#02)  compares  with  Robey’s  equation  (18)  for  the 


radiation  impedance  of  a  single  uniformly  pulsating  cylinder  in  an 
infinitely  long,  rigid  cylindrical  baffle  if  the  variables  of  integration 
are  changed  and  it  is  allowed  that  Robey’s  equation  contains  an  error  in 
sign#  (This  is  demonstrated  in  Appendix  II)# 

It  will  now  be  shown  by  uso  of  the  Dirac  Delta  Function  that 
equation  (3*02)  will  transform  into  the  well  known  expression  for  the 
radiation  impedance  of  an  infinitely  long,  uniformly  pulsating  cylinder# 
Equation  (3*02)  is  written  in  the  form  of  equation  (2#19)  by  trans¬ 
formation  of  variables: 

H#  Robey,  On  the  Radiation  Impedance  of  Finite  Cylinders,  J.  Ac  oust. 
Soc.  Am#,  27,  pp.  706*711,  1955 
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^oo 


<2>  = 


sin  (2  7 T 


Ho^fiT'ka) 


-oo 


^(2^  )yV„  (2) 


dr 


H,  (<Tka) 


(3.03) 


Consider  the  expression 

liM  1  ijr) 

2»y»*oo  2T//3  (^  )2 


(3.04) 


Expression  (3.04)  has  the  following  qualities: 

1#  It  is  zero  for  all  j°  except  J3  -  0. 

2.  At  f  -  0,  (3o04)  tends  to  infinity. 

3.  The  integral  of  (3.04)  from  J3  equals  minus  infinity  to 
^  equals  plus  infinity  is  equal  to  unity. 

The  above  is  tho  definition  of  a  Dirac  Delta  Function.  Therefore: 

.  2, 

1  S‘“ 

Lm  — *  - 

2  Up- 

(3.05) 

and 


sinf(£S^J  .  ^ 

7r  ^  v  cl  y 

2  7T-/3  (  r  )? 


lim 

i ->oo  ^  o  v/2  i  p  ^2 


oo 


lim  <z> 

2?Tj3^oo 


=  ^oc 


-oo 


k  ( JT) 


(erka) 

(2)  '  ~ 

H,  ({T“  ka) 


d  r 

J 


H 


(2) 


(ka) 


lm  <Z>  =  j£o0  -r-r 

■A— k  no  .. 


2lT|'i~*-  OO 


H, 


(ka) 


(3.07) 
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The  equation  (3.07)  which  is  the  radiation  impedance  per  unit  area 
of  an  infinitely  long,  uniformly  pulsating  cylinder  compares  with  the 
well  known  expression  for  such  radiation  impedance.®  This  concludes  the 
discussion  of  the  validity  of  equation  (2.27). 

£ 

L#  Beranek,  Acoustic  Measurements,  John  Wiley  and  Sons,  Inc#,  1949, 

p0  60 
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4.  Beam  Steering 


The  application  of  equation  (2# 27)  to 5 specific  array  will  new  be 
considered  in  order  to  study  the  effect  of  electrical  bean  steering  on  the 
radiation  impedance  of  the  array.  The  array  is  to  be  constructed  as  follows 


1.  There  aro  nine  active  cylinders  with  kp  -  -g-TT',  ka  and 


kb  =  0 

2.  The  cylinders  are  mounted  coaxially  in  an  infinitely  long, 
rigid  cylindrical  baffle  with  ka  -  Y 

3.  The  cylinders  pulsato  uniformly  over  their  surface  (symmetry 
in  azimuth). 

It  is  proposed  to  investigate  the  effect  of  electrical  beam  steering  on 
the  radiation  impedance  of  the  individual  cylinders.  Electrical  beam 
steering  will  be  accomplished  by  phase  leading  each  cylinder  a  uniform 
number  of  degrees  (  i )  with  respect  to  its  position  in  the  array. 

Equation  (2.27)  is  written  to  conform  to  the  above  boundary 
conditionso 


r 

02 


9 


<2>h  =  j 


0 


-/  j 


0 
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(4.01) 


The  result  of  the  evaluation  of  (4*01 )  by  numerical  integration 
using  the  trapezoidal  rule  and  A  f  =  — -  is  plotted  in  figures  2  and  3 
for  cylinders  number  one,  three,  and  five# 

For  such  an  array,  the  electrical  beam  steering  has  but  a  moderate 
effect  on  the  radiation  resistance  of  any  element.  Electrical  beam 
steering  affects  the  radiation  resistance  of  tho  end  elements  of  tho  array 
to  a  greater  degree  than  it  affects  that  of  the  center  elements. 

Electrical  beam  steering  has  a  more  drastic  effect  on  the  radiation 
reactance  of  such  an  array#  Radiation  reactance  in  this  array  varies  as 
much  as  ten  to  one  for  the  center  element  of  the  array# 
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APPENDIX  1 


The  Helmholtz  Wave  Equation  is  restated  in  cylindrical  coordinates: 

IM.  /A  JLl  yi  4  0  (AX. 01) 


9  r^~  ^  r  3  r  ^  £  e^  7  3  z2 


The  solution  to  the  above  equation  is  assumed  to  exist  in  a  variables- 
separable  form: 

0(r,e,z)  =  R(r)9(e)Z(z)  (A1.02) 

(Al#02)  is  substituted  into  (A1*01)  and  the  resulting  equation  is: 

1  d^R 
R  dr‘ 


/  1_  _dR_  ,1_  d2fr  ,  1  d2Z  ,2 

*2  Rr  dr  ftr2  d«2  Z  dz^ 


(A1.03 ) 


1  d^Z 

It  can  now  be  seen  that  —  — *  is  equal  to  a  quantity  which  is  independent 

Z  dz^ 

of  z;  therefore  the  following  is  true: 


-i — _  _  (2  ■fr  ^  )2 


(A1.04) 


Z  dz2 

Where  Jp  is  an  eigenvalue  to  be  determined  by  the  boundary  conditions 
assigned  to  the  array*  The  solution  for  Z(z)  may  now  be  written# 


Z(z)  =  eJZ'K'S  z  -/by  e  z 

Continuing  with  the  solution,  it  can  be  seen  that 

to  a  quantity  that  is  independent  of  e;  therefore: 

-L  =  .  n2 

9  d9 


d29 


9  de 


(A1.05) 
is  equal 


(A1.06) 

Where  n  is  an  eigenvalue  to  be  determined  by  the  boundary  conditions 
assigned  to  the  array*  It  will  be  found  more  convenient  to  write  the 
solution  for  9(e)  as  follows: 

9(e)  "  cncos(ne)  ■/  ci^  sin(ne)  (Al*07) 
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The  assignment  of  arbitrary  eigenvalues  to  the  solutions  for  Z(z)  and  9(e) 
determine  the  form  of  solution  which  (r)  must  assume# 

d\fy  1  <jjk2  _  (21 r'S  )2  ]  n2  j,  _0  (A1.08) 

dr  r  dr  r 

In  the  following  discussion,  the  quantity  k2  -  (2ft^  )2  will  be  defined  as 
equal  tO/^~^«  The  quantity  /^  must  reduce  to  k  as  (2^^  )2  approaches  zero 
therefore : 

=  ^k2  -  (2 ~7r'f  )2'  (2 TTf)z  <  k2 

-j  J(27r^  )2  -k2  (27?J)2  >  k2  (A1.09) 

To  satisfy  the  Sommerfeld  Condition  of  Radiation,  the  solution  for  Rn^  (r) 
must  be : 

Rn> (r)  =  Hn(2)  <yur)  (A1.10) 

The  expression  for  Rn  ^  (r)  as  r  approaches  infinity  is: 

R  (r)  r  1  2  e_j(/  r  (Al.ll ) 

This  expression  coupled  with  the  time  dependence  eS14^  represents  an 


out  g  o  ing  wave  # 

The  wave  equation  has  been  solved  consistent  with  the  Sommerfeld 
Condition  of  Radiation#  It  is  now  possible  to  state  a  eigenfunction  which 
is  a  solution  to  the  wave  equation;  and  if  the  boundary  conditions  imposed 
on  the  array  are  consistent,  a  sum  of  the  eigenfunctions  should  satisfy 
the  additional  boundary  conditions# 

The  eigenfunction  is  stated: 


a*  <*•••*>■*»*  :“fcKJ2"r2Hn(2Vr) 


(A1.12) 
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APPENDIX  II 


To  substantiate  the  statement  that  equation  (18)  of  Robey's  paper  is  in¬ 
correct,  the  development  of  equation  (18)  from  equation  (16)  will  be  inves¬ 
tigated  as  yrell  as  the  validity  of  equation  (16).  Further  it  will  be  shown 
that  equations  (19)  and  (20)  of  Robey's  paper  are  correct. 


zqq  = 


,oo 


H0(2)  (Ka)  sin2(jL) 
K  H-^^Ka) 


dj 


(16) 


Where  K  =  (k2-f  2)« 
Rewrite  equation  (16)  as  follows 

7^p° 

4a  feu) 


4^L 


sm 


'll  L) 


H, 


(2) 


(Ka) 


(ff  K  Hp(2)  (Ka) 


df  (A2.01) 


-oo 


Where  <  z^q>  =  *th e  average  radiation  impedance  per  unit  area 
Consider  the  following  expression: 


lira  BAn-2ILLi 
L-^oo  ff  L(  j°  )2 


(A2.02) 


The  following  facts  are  true: 

1.  (A2.02)  approaches  infinity  as  Jp  approaches  0. 

2.  (A2.02)  is  zero  for  all  other  f  . 

3.  The  integral  of  (A2.02)  from  J5  =  -oo  to  6P  a  -/oo  is  equal  to  unity. 
Therefore  (A2.02)  is  a  Dirac  Delta  Function,  and  (A2.01)  may  be  written 

as  follows: 


lim  <zqq>  * 


L~>oo 


H0(Z)(Ka) 
K  II1(2)(Ka) 


(A2.03) 


-oo 
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Thus  the  radiation  impedance  per  unit  area  of  tho  uniformly 
pulsating  cylinder  as  derived  by  Robey  reduces  to  the  following  as  the 
length  of  the  cylinder  approaches  infinity: 


iim  Czqq>  =  -jfo° 
L-*  oo 


Hq(2)  (ka) 

Hj<2)  (ka) 


(A2.04) 


This  does  not  agree  with  the  well  known  result  for  an  infinite 

IIQ(2)(ka)  . 


cylinder*  Further  the  angle  associated  with  "Q"/ '  is  always  in  the 

III  (^) 

fourth  quadrant;  thus  the  radiation  resistance  given  by  (A2#04)  would 
always  be  negative# 

It  oan  be  assumed  that  Robey’s  equation  (16)  is  correot  except  for 
sign,  and  it  will  be  shown  that  Robey’s  equation  (18)  has  an  error  in 
sign  but  that  equations  (19)  and  (20)  are  oorrect#  Equation  (16)  is 
restated  with  the  correct  sign: 


zqq 


sin^  (J°L)  Hc(  ^  (Ka) 
</>2  ™i  <2>  (Ka) 


(A2.05) 


zqq  =  d8a^ 


J 


sin2  (f  L)  Hq(2)(  f  k2-y2’  a) 

(1  )  2  |k2-f  2'  Ilj^H  |k2  -  fS  a) 


4 


k 


sin2(  J°  L) 
</)2 


(-j)  K0(  \ 

J2-k2'  a) 

(-j)\|/2-l 

^  Kl(  ^/2-k2'a) 

dy 


(A2.07) 
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To  reduce  (A2.07)  to  tho  form  of  Robey's  equation  (18 )  consider  the 


following  integrals: 
k 

*1  = 


Sin2(y  L)  Ho(2)(  I**-  X 


(f  ) 


Hl{2,(  [7-  J>2'a) 


dj> 


Let  'f  =  k  sin  T?' 
then  when  =  0,  "tf"  s  0 

Jn  k,  -XT  -  nf 
2 


and  |k2  ^  Cos^ 

dy  a  k  cos-^  dz*' 

Therefore 


'2 


sin^(kL  sinz?)  H0^^(ka  cost?") 


(siniy)' 


(ka  cos  t?) 


d  & 


Consider 


*2= 


oo 


sin2(J>  L) 


U 

k 


(jr 


K0((y2  -k2  a) 


2-k2  a) 


Let  y  :  k  cosh/7 

then  when  If  :  k,  «  0 
y  =  oo,  oo 

and  Uz  -k^  =  k  sinh^ 

df  =  k  sinh^  d<y 


df 


(A2.07) 


(A2.08) 


(A2.09) 
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Therefore 


o© 


X2  =  ~*2 


sin2(kL  cosh  )  Kp  (ka  sinh  ) 

(cosh-/')2  Kx  (ka  sinh  ) 


d  f 


(A2.10) 


Equation  (A2.05)  may  now  be  written  as 

df 
'2 


zqq  =  J 


8  a  (j-) 


r 

j 


o 

/-  o 


8a  'Pncu 


) 


2  (2) 
sin  (kL  sini^)  1IQ  y  (ka  cos  ) 

(sinx*)2  oos^) 


sin2(kL  coshy/)  Kp(ka  sinh^'  ) 
(cosh/'  )2  K-^ka  sinh/-) 


■/ 


d'f  (A2.ll ) 


Thus  from  (A2.ll),  it  can  be  seen  that  the  terms  in  Robey's  equation 
(18)  should  have  identical  signs;  and  from  the  initial  argument  about 
equation  (16),* the  terms  should  be  positive. 

Robey's  expressions  for  $  and &  as  they  appear  in  equations  (19)  and 
(20)  will  now  be  derived  from  (A2.ll).  The  C/  torn  is  self  explanatory  and 


correct.  Examine  the  expression 

Hn(2)(ka  oos  "O'  )  a  Hq(2)  (x) 


II- 


•j_  ^2  ^  (ka  cos  ~iy  ) 


H- 


(2) 


where  x  -  ka  cos 


(x) 


Jp(x)  -  jYo(x) 
J].(x)  -  JY1(x) 


Ji(x)  4  JYiU) 


JqU)  /  Y]_(x) 

JQ(x)J1(x)  7^Y0(x)Y1(x)  ^  J0(x)Yi(x)  -  J1(x)Y0(x) 

J2(x)  /Y2(x)  J1(x)/Y1(x) 
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J0(ka  Gos^)Ji(ka  cos  /Y0(ka  cos^)  Y\  (ka  cos  <^) 


J 


ka  cos 


2 


Thus  it  can  be  seen  that  Robey1  s  coefficients  Jc  ,  $  ,  and  £  in 
equations  (19)  and  (20)  are  correct  as  stated  but  are  not  consistent  with 
equation  (18)o 

The  above  arguments  are  believed  to  be  sufficient  to  show  Robey1 s 
equation  (18)  is  not  correct  as  stated,  but  would  be  correct  in  the  form 
of  equation  (A2.ll)*  Equation  (A2*ll),  if  stated  in  terms  of  average 
radiation  impedance  and  with  a  change  of  variable,  is  equivalent  to 
equation  (3#02)# 

7 

Go  N.  Watson,  A  Treatise  on  the  Theory  of  Bessel  Functions,  The  MacLiil 
Company,  1945,  p  77 
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(Figure  J.  ) 
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